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Abstract 

We show that the next-to-leading-order renormalizat ion-group- improved asymptotically- 
free BFKL Pomeron provides a good fit to HERA data on virtual photoproduction at small x 
and large Q 2 . The leading discrete Pomeron pole reproduces qualitatively the Q 2 dependence 
of the HERA data for x ~ 10~ 3 , and a fit using the three leading discrete singularities 
reproduces quantitatively the Q 2 and x dependence of the HERA data for x < 10~ 2 . This 
fit fixes the phase for all the BFKL wavefunctions at a chosen infrared scale. 
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The nature of the QCD Pomeron continues to perplex and intrigue both experimentalists 
and theorists [lj. Study of the Pomeron has been one of the most interesting aspects of 
the HERA experimental programme, with the discovery of a 'hard' Pomeron in virtual 
photoproduction [2] whose relation to the 'soft' Pomeron that is familiar from traditional 
hadronic reactions [3] is still the subject of theoretical speculation. In the near future, the 
LHC will provide possibilities to test theoretical approaches that have been honed with 
HERA data and may provide novel opportunities to study new Pomeron physics. 

Most of the HERA data on deep-inelastic structure functions are described well by 
the asymptotically-free renormalizat ion-group evolution expressed in the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) [4j equations. On the other hand, it has been suggested 
that a more appropriate framework for describing data at very low x is the diffusion in trans- 
verse momentum incarnated in the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation [5]. 
There has been considerable discussion of HERA data at low x in the context of unified 
BFKL and DGLAP equations [6]. However, it has not yet been established whether the pure 
BFKL Pomeron provides an accurate description of data in the kinematic range accessible 
to HERA. 

The BFKL equation with fixed strong coupling yields a leading Pomeron singularity that 
is a Regge cut, not a pole. Fixing the QCD coupling may be plausible in a suitable infra-red 
limit, but the coupling runs significantly in the ranges of virtuality and transverse parton 
momenta explored in inclusive measurements of structure functions at HERA. Over 20 years 
ago, it was shown [7] that within the BFKL formalism, the leading Pomeron singularity is 
a discrete Regge pole if the strong QCD coupling is treated correctly as asymptotically free 
and the infrared behaviour is encoded in a fixed phase. This leading singularity is followed 
by an infinite series of lower-lying poles. 

The next-to-leading-order (NLO) QCD correction to the leading-order (LO) asymptotical- 
ly-free BFKL equation is known [Sj, and it has been shown how to re-sum higher-order 
corrections so as to tame the NLO corrections [9j. The asymptotically- free BFKL equation 
described in [7J describes both the x- dependence of the un-integrated gluon distribution and 
the dependence on the transverse momentum, k, of the gluon, and hence also the Q 2 depen- 
dence of the structure function at low x. The only unknown quantity is the proton impact 
factor, which describes the couplings of the proton to the Pomeron trajectories. One 

benefit of the discrete Pomeron approach is that a simple expression for the un-integrated 
gluon density can be obtained, in terms of a small number of parameters, by expanding 
the proton impact factor in terms of the discrete set of solutions of the asymptotically-free 
BFKL equation. We recall that the BFKL Pomeron may be expressed as an integral that 
includes DGLAP as a saddle-point approximation valid in the double limit ln(l/x) ^> 1 and 
a s (Q 2 )\n(l/x) <C 1 [10J. However, this DGLAP approximation to the BFKL integral is no 
longer valid when the second condition is not satisfied, as in the case of low-x HERA data, 
where the discrete series of BFKL Pomeron Regge poles is a better systematic approximation 
scheme. 

We are unaware of any overall fit to HERA data made using the re-summed NLO 
asymptotically-free BFKL Pomeron. We perform such a fit in this paper, and show that it 
describes the inclusive virtual photoproduction HERA data both qualitatively and quantita- 
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tively. We show first that the leading BFKL Pomeron pole provides a successful qualitative 
description of HERA data on inclusive virtual photoproduction at small x ~ 10~ 3 , over a 
large range of Q 2 . This fit improves if lower- lying BFKL Pomeron poles are included, and 
we show that the asymptotically-free BFKL approach provides an excellent quantitative fit 
to all the inclusive HERA data at x < 10~ 2 , if the three leading BFKL Pomeron Regge poles 
are included. As well as the residues of the three BFKL Pomeron poles at zero momentum 
transfer, t, the BFKL fit has an additional free parameter corresponding to the value of the 
phase of the BFKL wave function that is assumed to be fixed by infra-red dynamics at a 
momentum k ~ 0.3 GeV0. 

We consider first the BFKL analysis of a zero-momentum-transfer process, at fixed strong 
coupling, a s . In this case, the eigenfunctions of the BFKL kernel are representations of the 
two-dimensional conformal group in the space of the transverse coordinates of the gluons, p. 
We include the BFKL characteristic function up to NLO [8], and use the re-summation of 
Scheme 3 proposed by Salam [9], which moderates the correction to the leading intercept as 
well as preserving the sign of the curvature of the characteristic function near the intercept, 
up to large values of a s . 

Considering only the leading conformal spin, the eigenfunctions may be written in mo- 
mentum space as 

w>) = S§, ID 



with 

Uk) =(k 2 ) l \ (2) 
where the eigenvalue uj is the solution to the equation 

uj = x(a«,f) = a s (l-Aa e )35^- + a s B + ^+iv^+afcfi(v). (3) 



Here 

- = 9± 

Qt s — Otsi 
TV 



Xb(*)=2(V>(l)-3te Mz)]) 



and 




where Ca = 3 and n/ is the number of active flavours at momentum k. Xi{ u ) is the 
NLO characteristic function given in [8], omitting the conformal symmetry- violating part 
associated with the running of the coupling (which is subtracted so that the 0(a 2 s ) terms on 
the RHS of (j3J) are not double-counted: see [9]). The implicit equation (j3J) for uj is readily 
solved using an appropriate combination of Newton's method and iteration. 

1 The precise value of ko is not an essential parameter. 
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Turning now to the case of running coupling, it was shown in [7] that the frequency v 
of the oscillations acquires a dependence on k, such that for a fixed eigenvalue u, u u (k) = 
X^ 1 (uj,a s (k)) is the solution to 

u =xMk),u UJ (k)). (4) 

This leads immediately to a critical value of the transverse momentum, fc crit , such that 

u = x(a s (^crit),0). (5) 

Provided x"( a s(k crit ), 0) is negatived, the value of u u (k) becomes imaginary for k > k CT it 
and the eigenfunction decreases exponentially as A; — ► oo. It is in order to ensure that 
x" { a s{k cr it), 0) remains negative that we have opted for Scheme 3 of the re-summation pro- 
cedure described in [9]. 



For k ~ fc C rit, the BFKL equation may be approximated as 



d 2 



+ 



fa X(«s(^crit),0) 



In 



k 2 



with 



dpn(*7*lit)] 2 2vr X "(a s (A; crit ),0) \k, 

11C A 2 



2 

crit . 



fM = 0, 



(6) 



0o 



3 3 

We recognize this as Airy's equation with argument proportional to ln(k 2 / k^ rit ) . Away from 
^crit; provided the running of the coupling is not too fast, so that 

« -w. 

the BFKL equation may be approximated semi-classically by 

d 



which has the solutions 
where 



dln(k 2 ) 



Vuj(k) 



Uk) = 0, 



(7) 



fc crit dk! 
k' 



In all regions, the solutions decrease as k — > oo, and are well approximated by 



Uk) = V3^ ip w (k)Ki (tp w (h)) (k > k ciit ), 



whereas 



f u {k) = y<p u (k) Ji(tp u (k)) + J_i(y u (fc)) (k < fccrft), 



(9) 

(10) 
(11) 



We use the notations x"(a s ,f) = d x(a s , and x(a s ,^) = dx(a s , v)/da s . 



4 



where we have expressed the appropriate Airy function in terms of the modified Bessel 
function of the second kind, Ki, and Bessel functions of the first kind, J ± i. Away from 
k ~ k crit where <p u becomes large, these Bessel function solutions approximate the solution 
to the semi-classical equation (jTJ). 

It is important to note that the matching of the solutions at k — k crit determines the 
phase of the oscillations in the region where k < k CTit , for a given value of uj. Following [7J, we 
encode the unknown infrared behaviour of QCD by assuming that it leads to a fixed phase, 
77, at some low value of the transverse momentum, ko, which we take for definiteness to be 
0.3 GeV@. More precisely, the infrared condition is given by 



/■ fc crit dk' ( 1\ 

<pM = 2^ -yMky = (n-~)ir + r ] , (12) 



rfccrit dk' 
lk 

and means that, just above k = ko, the wavefunction behaves like 



Uk) ~ sin (k' - kl) - r\ . (13) 

Once the phase condition (fT2l) is imposed, only a discrete set of values of the eigenvalue uj are 
allowed simultaneously by the infrared phase condition and the phase condition imposed by 
the matching, giving rise to a description of the QCD Pomeron as a discrete set of isolated 
poles, as opposed to the cut found if the running of the strong coupling is neglected. 

In order to express the low- 2 structure function of the proton, F 2 (x,Q 2 ), in terms of 
these eigenfunctions, the eigenfunctions themselves must be convoluted with the impact 
factor $ p (k), that describes how the proton couples to these trajectories at zero momentum 
transfer. In the case of the un-integrated gluon density xg(x, k), we have 

r dk' 

xg(x,k) =Y.j jr%(k')x-^k 2 f: n (k')U n (k), (14) 
and the un-integrated gluon density is related to the structure function by 

r Q dk 

F 2 (x,Q 2 ) = / — <S> Dls {Q,k)xg(x,k), (15) 
Jo k 

where the impact factor, $dis ; that describes the coupling of the virtual photon to the 
trajectories is given by (see [1]) 

* ^ m ^2 ^ 2 f 1 , , 1 " 2p(l - P) - 2r(l - r) + 12p(l - p)r(l - r) 

(16) 

The proton impact factor, <& p (k) is unknown a priori and has to be fit to data. Since the 
eigenfunctions f Un (k) form an orthonormal set, we can expand the impact factor as a series 
in these eigenfunctions with a discrete set of coefficients, a n : 

%{k) =Y / a n k 2 f uln (k), (17) 



3 A change in this value of the infrared momentum scale can be compensated by a change in the phase, 
77, so that the infrared behaviour of QCD is in fact encoded using a single parameter. 
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and exploit the orthogonality properties to write 



xg(x,k) = ^ 



a n x n k ftj n (k). 



(18) 



n 



A model for could be used to estimate the coefficients a n , which could be also con- 

strained using other HERA data, e.g., on the diffractive production of vector mesons. 

At sufficiently small x, we expect this sum to be dominated by the first few poles. The 
contribution from the remaining poles could be approximated by assuming that the effect 
of fixing the phase at k on the allowed values of us is negligible for uj < 0.1, and that the 
discrete set of eigenfunctions may be replaced by a continuum. In this case, one simply 
adds to the expression ( fT8|) for the un-integrated gluon density the following integral that 
represents the contribution from such a continuum: 



where bpiy) is a function that encodes the coupling of the proton to all the remaining eigen- 
functions, and is chosen to be real so that the wavefunctions in this continuum also respect 
the imposed infrared condition ([TBI . In order to implement such a programme one would 
need additional parameters to characterize the arbitrary function b p {v). However, in the re- 
gion of x and Q 2 considered, we find an excellent fit without making use of such a continuum 
and the associated extra parameters, and hence do not consider it further. Nevertheless, it 
should be emphasized that at Q 2 that is sufficiently large for the DGLAP analysis to become 
valid, the double-leading-logarithm DGLAP behaviour would be embedded mainly within 
this continuum contribution. 

We have determined numerically the eigenfunctions of the leading four poles of the NLO 
asymptotically-free BFKL Pomeron. We limited ourselves to the first four poles because 
their u values are in the same range as the observed rate of rise, A, of the F% measurements. 
This is determined by fitting the measured F2 to x~ x at fixed Q 2 and is closely related to the 
logarithmic derivative <ilog(i 7 2)/(ilog(l/x). The values of A determined phenomenologically 
by experiment vary between A « 0.1 for Q 2 < 0.6 GeV 2 and A « 0.33 for Q 2 > 60 GeV 2 [UJ. 
The leading eigenvalue, u±, depends on the infrared phase, 77, varying between u\ = 0.235 at 
77 = and u)\ = 0.315 at rj = tt/2, The sub-leading eigenvalues are smaller, the fourth one, 
co>4, being ps 0.10. 

We determined simultaneously the best-fit value of the infrared phase 77 and the coeffi- 
cients, a n . The fit was performed in the low-x region, x < 0.01 and for Q 2 > 4 GeV 2 , so 
as to avoid saturation effects. The saturation scale at HERA was found to be Q 2 ~ 0.5 
GeV 2 [12], implying that saturation effects should fall below the measurement precision for 
Q 2 > 4 GeV 2 [13]. The best fit is obtained for rj = — 0.217T, with the values for the first four 
eigenvalues and their corresponding fc cr it, given in Tabled]; the corresponding eigenfunctions 
(normalized in the domain k > ko) are shown in Fig. [TJ We see that the eigenvalues indeed 
decrease as n increases, so that for sufficiently small x the leading trajectories should be 
sufficient to describe the data over any fixed range in k. We note also that the eigenvalues 
approach each other as n increases. 




(19) 
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n 


UJ 


fccrit (GeV) 


1 


0.26 


5.9 


2 


0.17 


330 


3 


0.13 


2.8 x 10 4 


4 


0.10 


2.6 x 10 6 



Table 1: The eigenvalues and values of k cr it for the 4 leading eigenf unctions of the 
asymptotically-free BFKL Pomeron, for n = — 0.217T at ko = 0.3 GeV. 
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Figure 1: The first four eigenf unctions of the NLO BFKL kernel with running coupling and 
infrared phase rj = — 0.21-7T at ko = 0.3 GeV. The arrows indicate the values of k crit . 



Number of poles 


X 2 /N df 


a\ 


a 2 


a 3 




1 


3624/101 


0.035 








2 


264/100 


0.029 


-0.028 






3 


91.4/99 


0.041 


0.055 


0.085 




4 


91.3/98 


0.042 


0.067 


0.11 


0.016 



Table 2: The qualities of fits using up to 4 poles, and the corresponding pole residues, as- 
suming r\ = — 0.217T at ko = 0.3 GeV. 
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Having fixed the value of 77, we investigate the number of eigenf unctions required for 
a good description of the HERA data for x < 10~ 2 . An overall 1-pole fit using only the 
leading eigenfunction has very poor quality: \ 2 l^df = 3624/101, though it does reproduce 
qualitatively the data for x ~ 10~ 3 , where it is more likely to dominate over the non-leading 
Pomeron poles. The quality of the overall fit improves significantly when the two first 
eigenf unctions are used: x 2 /Njf = 264/100, and the 3-pole fit is excellent: x 2 /N$ = 91.4/99. 
On the other hand, adding a fourth eigenfunction does not improve the fit any further: 
X 2 /N d j = 91.3/98. Since also the coefficient of the leading eigenfunction, ai, is almost the 
same in the 3- and 4-pole fits, in the following we consider only the fits with 3 or less 
eigenfunctions. 

Fig. [2] compares the results of the 1- and 3-pole fits with the measured values of F 2 . We 
see that the 3-pole fit indeed describes the data very well, corresponding to its very good 
X 2 ■ Fig. [2] also displays the 1-pole fit; despite its very high x 2 , it reproduces qualitatively 
the main features of the data, particularly for moderate Q 2 . We note that the coefficient, 
di, of the leading eigenfunction in the 1-pole fit is about 20% smaller than in the 3-pole 
fit. This indicates that the excellent agreement of the 3-pole fit with the data is due in 
part to cancellations between the different eigenfunctions. To illustrate the properties of the 
3-pole fit, we show in Fig. [3] the contributions to F 2 from the 3 eigenfunctions separately 
as functions of the momenta k 2 at several characteristic x values. Fig. [3] shows that, in 
the region of medium Q 2 values: 4 < Q 2 < 20 GeV 2 , the contribution of the leading 
eigenfunction coincides with the fitted F 2 curve, i.e., the contributions of the second and third 
eigenfunctions cancel each other. However, at larger Q 2 , especially above Q 2 > 100 GeV 2 , 
the fit has large cancellations between all three components, and the leading eigenfunction 
cannot fit the data by itself. 

Fig. H] compares the Q 2 dependence of the effective value of the exponent A determined 
from a phenomenological fit to the data, and as extracted from our fits @. In the case of the 
1-pole fit (dashed line), A is identical with the leading eigenvalue: uj = 0.26, and in the 2-pole 
fit (dotted line) the values of A become mostly smaller than the leading eigenvalue. However, 
in the 3-pole fit, whilst the A values are smaller than the leading eigenvalue for Q 2 < 20 
GeV 2 (solid line), they become larger at higher Q 2 , and are closer to the values extracted 
form a phenomenological fit to the data. The surprising fact that the sum of the contribution 
with small eigenvalues can give a larger rate of rise than the leading eigenvalue is due to 
the fact that A is closely connected to the logarithmic derivative, A ~ rflog(F 2 )/<ilog(l/2;). 
Owing to these cancellations, the logarithmic derivative can become larger than the largest 
eigenvalue. The fact that the 3-pole BFKL fit gives somewhat smaller values of A than 
the phenomenological fit for Q 2 ~ 20 to 70 GeV 2 is closely related to the fact that the 
lowest-x point at each of these values of Q 2 lies slightly above the 3-pole fit, as seen in the 
corresponding panels of Fig. [2J However, the minor discrepancies for these few points does 
not spoil the quality of the overall fit, which is a better measure of its validity than the A 
plot shown in Fig. HI 

In summary: we obtain a very good description of the HERA low-x data in a large range 



4 In Fig. [4] we have also included recent data from Zeus [14] which are fully consistent with previous data, 
and therefore have not been used in the fit. 
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Figure 2: Comparison of the HERA F2 data with the 1- and 3 -pole fits, shown as dashed and 
solid blue lines, respectively. 
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Figure 3: The contributions to F 2 of the three eigenf unctions of the 3-pole fit. 
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Figure 4: The rate of rise X, defined by F 2 oc (l/x) x at fixed Q 2 , as determined in the three 
fits and in a direct phenomenological fit to the data. 
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of Q 2 , from 4 < Q 2 < 650 GeV 2 , using just three eigenfunctions and adjust 4 free constants: 
the phase 77 and the coefficients 01,2,3 □• Important roles are played in the fit not only by 
the leading eigenfunction, but also by the pattern of cancellations between the sub-leading 
trajectories, which is very sensitive to the parameter 77. For this reason, the quality of the fit 
is also very sensitive to the value of 77. Thus, for 77 = — 0.37T the x 2 grows to 142 for the 3 pole 
fit (instead 91.4 at the minimum), and at the extreme values of 77 = — n/2 and 77 = the \ 2 
values are 430 and 680, respectively. Consequently, the data determine the infrared phase 
quite precisely (within the theoretical framework described above) 0: 77 = —0.21 ± 0.027T In 
turn, the leading eigenvalue is also precisely determined: u = 0.26 ±0.01. The relatively low 
value of this eigenvalue is responsible for the fact that, although the A plot is reproduced by 
the 3-pole fit only in a qualitative way, we obtain a very good overall fit to the data. 

To our knowledge, this is the first time that the discrete asymptotically-free BFKL 
Pomeron has been shown to fit the HERA data at low x and high Q 2 . As such, we be- 
lieve that it is also the first time that a parametrization of the Pomeron derived from first 
principles in QCD has been confronted successfully with experimental data. A natural next 
step would be to extend this comparison to include other low-x HERA data, including those 
on the diffractive production of vector mesons, etc. One could also envisage the development 
of a BFKL Pomeron calculus and its deployment to make predictions for both inclusive and 
exclusive phenomena at the LHC. 
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